Scaling of singular structures in extensional 
flow of dilute polymer solutions 



Paul Becherer^'*, Alexander N. Morozov^'^ and 
Wim van Saarloos ^ 

^Instituut-Lorentz for Theoretical Physics, Universiteit Leiden, Postbus 9506, 
NL-2300 RA Leiden, The Netherlands 

^School of Physics, University of Edinburgh, JCMB, King's Buildings, Mayfield 
Road, Edinburgh EH 9 3JZ, United Kingdom 



Abstract 

Recently singular solutions have been discovered in purely elongational flows of 
visco-elastic fluids. We surmise that these solutions are the mathematical struc- 
tures underlying the so-called birefringent strands seen experimentally. In order to 
facilitate future experimental studies of these we derive a number of asymptotic 
results for the scaling of the width and extension of the near-singular structures in 
the FENE-P model for polymers of flnite extensibility. 
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1 Introduction 



Recently, there has been renewed interest in the properties of extensional flows 
of dilute polymer solutions, in particular in a class of flows known as internal 
stagnation point flows, such as the four roll mill flow and the cross-slot or 
cross channel flows shown in Fig. [H Theoretical modelling of such flows has 
proven to be particularly challenging, and despite extensive experimental and 
theoretical investigation, some aspects of these flows remain to be elucidated. 
For example, Arratia et al. [Ij have recently found a bifurcation to asymmetric 
stationary flow patterns in a cross-channel flow of a dilute polymer solution, 
followed by a secondary instability to time-dependent asymmetric flows. In 
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spite of its simple appearance, this instability is not yet fully understood. 
Conversely, instabilities have been found numerically that have no clear ex- 
perimental counterpart (for example, Harris and Rallison PfS] and Xi and 
Graham [Ij). 

In traditional investigations of elongational flows, it has often been assumed, 
either implicitly or explicitly, that the base flow solution is the classical solu- 
tion [5j of the Oldroyd-B continuum equation in which the stresses are constant 
in space. However, Rallison and Hinch ^ already mention singular solutions 
of the UCM equations that are strongly peaked at the centre line along the 
outflow direction, and Renardy [7] recently pointed out that these singular so- 
lutions are also relevant for Wi < 1/2, where they do not actually diverge, but 
are still singular. Although such solutions are not easy to probe numerically, 
recent work by Thomases and Shelley |8] shows that these solutions do emerge 
spontaneously in high accuracy numerical simulations of a model problem for 
elongational flows. One reason for this may be that, as we will discuss, the 
constant stress solutions typically do not satisfy the physical boundary condi- 
tions at the edges of the flow region. A natural question that therefore emerges 
is whether these singular solutions are the mathematical counterpart of the 
so-called birefringent strands that have been found experimentally p|10llll] 




and numerically [T2|13lll4j . Such strands are thin regions of highly extended 

polymers along the central outgoing streamlines. A strand can modify the 
flow, leading to a "dip" in the observed velocity profile of the exit flow [l5fT6] . 
and it seems reasonable to assume that they affect any instabilities that might 
occur. 

The singular solutions mentioned by Renardy [7] arise in models like the UCM 
model which essentially assume that the polymers are infinitely extensible. 
A second important question that emerges from these studies is therefore 
how these singular solutions are modified when we consider finitely extensible 
polymers, as in the FENE-P model [T7]. The first step in this direction was 
made by Renardy [7j, who showed that for the Giesekus model, which limits 
the growth of elongational stresses, stress profiles remain smooth while stress 
gradients can diverge. 

The aim of this article therefore is to study the scaling of the main characteris- 
tics (e.g., the width) of these singular solutions strands in a simple way, and to 
analyse the modifications due to the finite extensibility. Our analysis is done 
for the ideal case of purely elongational fiow, which is a good approximation 
near the stagnation point. By confining the analysis to this simple case we are 
able to derive a number of explicit asymptotic scaling results which we hope 
will facilitate making the connection between the (almost) singular solutions 
and the birefringent strands. 

Our results are consistent with results for steady fiow by Renardy [7j, Thomases 
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(a) (b) 
Fig. 1. (a) The four-roll mill (b) Cross-channel flow. 

and Shelley [8] and Xi and Graham [4J. However, since our approach is essen- 
tially one-dimensional, we can obtain a better numerical resolution, and it 
may be hoped that some calculations that are not (yet) feasible in a full two- 
dimensional approach, such as eigenvalue calculations for stability analysis, 
may be performed within the framework we present here. It should be men- 
tioned, however, that the current experimental resolution |T8] is substantially 
worse than the numerical resolution, even for full 2-D numerics. 



The layout of this paper is the following. In Section 2 we summarize the 
equations for extensional flow of a UCM fluid, and the structure of the singular 
solutions discovered recently. After analysing extensional flow of a FENE-P 
fluid in Section 3, we derive in Section 4 various asymptotic results for this 
case. We end the paper with a brief discussion of the robustness of our results. 



2 Purely extensional flow of a UCM fluid 



Putting the stagnation point at {x,y) = (0,0), a purely extensional flow field 
is given by 

V = {v^,Vy) = e{x,-y), (1) 

where e is the elongation rate. This flow field satisfies incompressibility, V -v = 
0. In a nondimensionalized formulation of the UCM model, the flow field 
becomes 

v=ix,-y), (2) 
and the constitutive equation for steady flow is [5] 

T + Wi [{v ■ V)T - {Vvf ■ T - T ■ (Vv)] =Vv + {Vvf. (3) 

Here T is the stress tensor, v the velocity, and Wi the Weissenberg number, 
which for purely extensional flow we define as 

Wi = eA, (4) 
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with A the relaxation time of the polymers. The momentum conservation 
equation for creeping flow is 

V ■ T - Vp = 0. (5) 



Solutions for steady flow are found by inserting the pure extensional flow into 
the constitutive equation: 

T.. + Wi X— ^-1/— ^-2T,.J =2, 
\ ox ay J 

T,,^,,d/^^y?^\.0, (6) 



(9T dT 

yy y. 

dx dy 



T,, + Wi ( - y-^ + 2T, J = -2. 



If we assume a spatially uniform stress field, these equations reduce to 

2 —2 

T... = ; T^y = 0; T„„ = . (7) 

1 - 2Wi ^ 1 + 2Wi ^ ' 



It is well-known that for Wi > 1/2, this solution diverges and becomes un- 
physical [5J. If we no longer require that the stresses be constant in space, 
the stress fields lose smoothness even below Wi = 1/2, as was pointed out by 
Renardy [7|. Around and above this value of Wi, there are large stress gradi- 
ents around y = 0, reminiscent of a birefringent strand. We shall derive these 
solutions, essentially following Renardy's presentation. 

We choose as flow domain the strip {x,y) G M x [—1, 1]. On this domain we 
assume pure extensional flow, see Fig. [21 We assume that the stresses depend 
only on y: 

Txx = T^x{y); T^y = Txyiy); Tyy = Tyy{y) = p{y), (8) 

such that momentum conservation is obeyed. We impose boundary conditions 
for the normal stresses on the "inflow" boundaries P^ : 

T,..(±l)=e; Tyy{±l)=r]. (9) 



The equations reduce to 

T..(y) + Wi(-yT^,(l/)-2T,,(y)) =2, 

T,y{y) + Wi{-yZy{y)) =0, (lo) 

Tyy{y)+Wi{-yT;yiy) + 2Tyyiy)) =-2. 
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Fig. 2. Uniform extensional flow on an infinite strip. 
The solution can be given in closed forn{l|: 



2 




l/Wi-2 



T 



XX 



T 



xy — 



l-2Wi 


-2 



l-2Wi 



2 




l/Wi+2 



(11) 



yy 



1 + 2Wi 



+ [V + 



l + 2Wi 



The first part of the xx and yy stress components is the same as for the uniform 
solution ([7]). The second part has to be added to make the solution consistent 
with the imposed boundary conditions. Due to the fractional exponents in this 
part, none of these solutions is smooth, (except when 1/Wi is an integer). 
This becomes important when approximating these functions using spectral 
methods. 

As Eq. (fTTj) clearly shows, the classical constant stress solution ([7]) only exists 
for very special boundary conditions. Indeed, the non-smooth terms in Eq. (ITT!) 
are present at every Weissenberg number unless one chooses ^ = 2/(1 — 2 Wi) 
and 7] = —2/(1 + 2 Wi). As Wi approaches 1/2, these terms create a narrow 
region of large extensional stress, qualitatively similar to a birefringent strand, 
as was already observed by Rallison and Hinch [6]. 

More interestingly, the purely elongational flow, Eq. ([T]) does not support any 
boundary value for the shear stress T^y other than T^y = at the bound- 
ary [T9|20] . Forcing T^y 7^ at the inflow would inevitably modify the velocity 
proflle, Eq. ([1]): in addition to the purely extensional flow fleld, it would ac- 
quire a shear component. The importance of this modiflcation will be discussed 
in Section 5 where we address the relevance of our results to the stability of 
experimental realizations of stagnation-point flows [I]. 

Note also that even in the range where the uniform solution clearly breaks 
down (Wi > 0.5), the non-uniform solution (fTT!) is well-deflned. However, the 



^ A similar solution was given in [7]. The UCM equations reported there, however, 
seem to be written for uniaxial extensional flow v = (1,— ^) rather than for 
planar extension v = (1,-1) as claimed in [Tj and as in the present paper. This 
explains the difference between Eq. pT]) here and Eq. (9) from [7]. 
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stress on the central outgoing streamline diverges, and for Wi > 1, the total 
elastic energy, which is proportional to the integral of trT, also diverges. 



3 Extensional flow in a FENE-P model 



The FENE-P model avoids the blow-up of the extensional stress by implement- 
ing a nonlinear force law for the polymer molecules [TT] . The solvent is usually 
treated explicitly. We shall ignore the solvent viscosity, thus formulating an 
extension of the UCM model rather than the Oldroyd-B model. 

We shall formulate the FENE-P model in terms of the conformation tensor 

A = {RR), (12) 

where R is the end-to-end vector connecting two beads of a single dumbbell. 
The brackets denote an ensemble average. 

In our nondimensionalization, the stress is given in terms of the conformation 
tensor as 

T- W ^ M (13) 

WiVl-trA/L2 1-2/^2^" ^ ^ 

Here, I is the unit tensor and is the maximal value of the trace of the con- 
formation tensor, that is, L is the maximal extension of the dumbbells, relative 
to their equilibrium extension. The 2 appears in the rightmost denominator 
because in two spatial dimensions this is tri. 

The constitutive equation for steady flow then is [8] 

f A I \ 



Wi 



V)A - {Vvy ■ A - A ■ (Vv) 



l-trA/L2 l-2/L2^ 



(14) 

Note that Eqs. and differ from the "classical" FENE-P model in 
that they are restricted to two dimensions and assume the other components 
of the conformation tensor {A^^, Ay^, A^z) to be zero. In the two-dimensional 
flow Eq. ([2]), the classical FENE-P model ^ would have A^z 7^ 0. However, 
this approximation bears no influence on the asymptotic result for the width 
of the birefringent strand found in Section 4. 

The momentum balance is nonlinear in the conformation tensor: 

V-A A-V(trA) , 

+ T^7^-rX77^ -Vp = 0. (15) 



Wi\l-trA/L2 L2(l -tr A/L2)2 
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Again, if we assume that the conformation tensor (and, hence, the stress ten- 
sor) depends only on this expression simphfies, and we find one equation 
for Ax;j. and Ayy, and an equation for A^y in terms of Ari;^ and Ayy-. 

1 (d^yy + ^yy ^^^All^] -Wi^ = 0. (16) 



and 



1 — trA/L^ y dy 1 — trA/L^ dy J dy 



1 (dA^y ^ A,y dilA/L^ \ ^ ^ ^^^^ 



l-tiA/L'^\dy l-trA/L2 dy 

There are a few things to note about this system. The first equation involves 
only Axx and Ayy. Suppose now that we can solve the constitutive equa- 
tion for Axx and Ayy] we can then always find a p{y) to satisfy the 
momentum balance. From the second equation we find that if A^y = is a 
solution of the constitutive equations (which in uniform extensional flow it is, 
as we shall see), then this always satisfies the momentum balance, similar to 
the UCM case. We can therefore restrict ourselves to finding a solution for the 
constitutive equations, given a uniform extensional flow. An analytical solution 
is no longer possible, and we will give numerical solutions. However, there is 
analytical information to be obtained from the equations, mostly asymptotics. 

Analogous to the UCM case, we assume that the extension depends only on 
?/, and we insert the pure extensional flow into the constitutive equation f[T^ : 



Wi{-yj^^M-'^AUy)) 
wi {-yKyiy) 

Wi[-yAyy{y) + 2Ayy{y) 



{AUy) + Ayy{y))lL^ l-2/L^ 

^xyiy) 



l-{AxM + Ayy{y))/L^' 



\l-{AUy) + Ayy{y))/L^ 1-2/LV 

(18) 

We shall use the same domain and boundary conditions as for the UCM case. 
For simplicity, we assume A^x = Ayy = 1 (the equilibrium values) at ?/ = ±1. It 
is clear that A^y = is indeed a solution of the equations. We can then ignore 
this component of the conformation tensor, and we may restrict ourselves to 
the diagonal components. 

For = 10, 100, 1000 we then find, for different Wi, the plots in Fig. [31 
These are clearly qualitatively similar to the birefringent strands that have 
been found experimentally and numerically Pll0flllll2lll3fl4j . 

The formulation in terms of a conformation tensor allows us to translate these 
results immediately to a birefringence profile: for the relative change in index 
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i^ = 10.0 Z.^= 100.0 1000.0 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

y y y 



Fig. 3. Extension for different L^, for Wi = 0.9,2.0,8.0, and Wi oo. Note 
the different scales on tlie vertical axis. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



y y y 

Fig. 4. Normal stress WiT^^ for different for Wi = 0.9,2.0,8.0 and Wi oo. 
Note the different scales on the vertical axis. The insets show a smaller range. 

of refraction n, we have the proportionality [21] 

^ oc ^{Ayy - A^^Y + AA%. (19) 

For the strongly stretched central region, the birefringence is approximately 
directly proportional to A^x- 

For the stresses, we can use Eq. (1131) . Since in our nondimensionalization, the 
physical stress is given by rj Wi T, with r] the viscosity, we plot WiT^j, rather 
than Txx in Figure |H 



4 Asymptotic results for a FENE-P fluid 

We present four types of asymptotic results. The first is the behaviour of the 
flanks of the strand. Here we shall find that we recover the UCM behaviour. 
Next, we look at the maximal "extension" of the strand — the value of 
-A^xx at 1/ = 0. It depends on Wi and L^. We then combine these results, and 
give an approximate expression for the width of the strand, which we define 
as the point where the UCM profile intersects A°^^. We also show that this 
gives practically the same results as another definition for the width, namely 
the inflection point of the Axx profile. Last, we look at the behaviour of the 
extension around ?/ = 0, and we show that even though the stresses at the 
centre stay finite, stress gradients may diverge for ?/ — 0. 
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L^=-f0.0 100.0 1000.0 




y y 



Fig. 5. Extension A^x, solid line, compared to "UCM extension" WiT^x + 1, dashed 
line, for different L^, for Wi = 0.9,2.0,8.0 and Wi — > oo. Note the different scales 
on both axes. 

4-1 Outer flanks 



For small extensions, the FENE dumbbells behave approximately as linear 
springs. We would therefore expect to recover UCM behaviour outside the 
centre of the strand, where the extension is relatively low. 

To compare the FENE-P profiles to the UCM results, note that from Eq. (fT3!) 
we find that for small extension at large 

T^—(A-1) for trA<L2 and > 1. (20) 
Wi 

For the xx component, this implies 

WiT,, (21) 

For the present boundary conditions, this means that outside the centre of the 
strand, we have 

Axx - WiT™ + 1 = (l - b|Vw.-) + , ^22) 

This approximation also works well for very high Wi, because the right hand 
side of Eq. (fT4|) is then negligible, until tr A comes very close to L^. Note 
that we compare the UCM stress with the FENE-P extension. For the small 
extension limit, the comparison of the stresses in both models also gives good 
results (in that case, the approximation fl2T]) holds well), but it breaks down 
in the large Wi limit. 

If we compare the profiles for the conformation tensor to the solution for the 
UCM fluid, as in Fig. 0, we see that for Weissenberg number above roughly 
1.0, the flanks of the stress profile are approximated by those of the UCM case, 
cut off in the centre by the value of Axx at y = 0. This observation helps us to 
obtain some asymptotic analytical results for the structure of the birefringent 
strand. 
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4-2 Maximal extension 



We can find the maximal extension by taking the limit of the equations (fT^ and 
( 1151) for y — s> 0. The advection terms then vanish, and we are left with the 

system 

-2 Wi A,^ = - ( , - -tAttt^] , (23a) 



^l-trA/L2 l-2/L2^ 

This system does not allow a simple analytical solution. However, we can make 
a fair approximation by assuming that A^x ^ Ayy and A^x ^ 1 (this requires 
Li^ ^ 1). We then have an autonomous equation for Axx'- 

This is easily solved and gives 

= - 2Wi) ■ 

This approximation can be shown to be self-consistent. For Wi > 2, clearly 
Axx ^ 1- Furthermore, since the flow is purely compressive in the y-direction, 
we expect Ayy < 1. In fact, we can find an asymptotic expression for Ayy at 
y = from Eq. ( ]23bl) . For large L^, the rightmost term is approximately 1, 



and assuming Ayy <^ Axx, the equation simplifies to 
Inserting the asymptotic expression for Axx, we find 



^^'^ ~ 4Wi' ^^'^^ 

For large L and Wi > 2, this implies that Ayy is of 0(1/ L"^) compared to Axx, 
and the above approximation is self-consistent. We conclude that for ^ 1 
and Wi > 2, Eq. (125|) is a reasonable approximation. 

This is confirmed by numerical solution of the system (!23l) . In Fig. [6] we show 
the maximum value of A^r as a function of Wi for different L^. 
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Fig. 6. Maximum of Axx as a function of Wi for different L^. Numerical solution, 
solid line, compared with the approximation (j25|) . dashed line. The arrow indicates 
the asymptotic value for Wi oo. Note the different scales on the vertical axis. 



4-3 Width of the birefringent strand 



Given these results, we can now proceed to derive an approximation for the 
width of the birefringent strand as a function of Wi and L^. We define the 
width as the location where the UCM approximation intersects the value of 
the "centre maximum" of the extension. This is essentially the width of the 
"plateau" in Fig. d 

We combine the results in eqs. flTTl) . and fl2^ to obtain for the intersection 
point Ho 

^(l-M—V 1^.^(1-^). (28, 
Solving for yo thus gives the strand half-width as 

Wi/(l-2Wi) 



yo 



1 9 / 1 \ 
+ 1 ^ 



2Wi V 2Wiy 



(29) 



For Wi — > 00 we find yo —>■ 1/L. An expansion around Wi = 00 gives 



1 2L2(1 -InL) - 1 1 ^ / 1 \ 



For ^ 1 this reduces to 

1 1 - InL 1 ^ / 1 \ 

The strand half-width in these approximations is shown in Fig. [7] for different 
values of L^. Note that the "exact" solution (l29l) bends up for lower Wi. This 
is unphysical. It shows up because both sides of Eq. (l28l) are approximations 
that break down for relatively small and Wi. 

The asymptotic value 1/L has already been found by Mackley and cowork- 
ers fTO\, in a more qualitative way. 
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Wi Wi Wi 

Fig. 7. Width of the birefringent strand as a function of Wi for different . "Exact" 
solution, Eq. (f29ll . solid line, compared with the approximation (fSTI) . dashed line. 
The arrow indicates the asymptotic value and the squares denote a numerical 
result where the width is defined by the inflection point of the stress profile. Note 
the different scales on the vertical axis. 

4-4 Central region 



The behaviour of the extension (and, hence, the normal stresses) as y ^ 0, 
is not immediately clear from the plots in Figures [3] and HI In particular, the 
question arises whether the stresses display singular behaviour for y — 0, as 
was found by Renardy for the Giesekus model [7]. 

It is clear from Eq. (fTSl) that ?/ = is a singular point of this system of 
differential equations: as ?/ — > 0, the higest-order derivative vanishes. The 
conformation tensor will therefore not be analytical in ?/ at y = 0. To determine 
the degree of this singularity, we follow the usual procedure for dealing with 
singular points [22], and we expand the solution around |/ = as 

oo 

A(|/) = ^a,/+'' + Ao, (32) 

r=0 

with (3 in general non-integer, and Aq the uniform solution of Eqs. (l23l) . We 
then seek an expression for the exponent (3 of the dominant (lowest-order) 
term by balancing the leading-order singular terms. For the small-extension 
limit, the equations actually reduce to the UCM equations, and the exponent 
is given by the exact solution Eq. (ITTl) . that is, 

(5 = — -2 for trA<L2. (33) 
Wi 



For moderate to high extension, we need to take into account the finite ex- 
tensibility of the polymers. To lowest nontrivial order in y, the expansion fl32l) 
becomes, in components. 



, (34) 
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Weissenberg number 



Fig. 8. The exponent /3 as a function of Wi. The dashed hnes are the asymptotic 
approximations: the curve is from Eq. (I33p . the two straight hnes are from Eq. (I38p . 



The dots were calculated from numerical solutions for a^^ and Uyy. 



We insert this into the constitutive equations flMl) for purely extensional flow, 
and expand the fractions in that equation to first order in a^x and ayy. We then 
have at first order, leaving out the constant terms (which form the equation 
for Ao), 



a. 



1 -trAo/L2 

a. 



+ + - 2 Wi a.. - /3 Wi a. 



"yy 



2^2 



tr Ao/L 

+ -^y^- + ""yyy^' +2Wiayy-(3 Wi ayy 



l-trAo/L2 (l-trAo/L2)2 



0, 
0. 



(35) 



We can simplify these equations considerably. We put Ayy 0, because ac- 
cording to (123bp it is of order l/L^ compared to the other terms and also 
compared to 1 — A°^^/L'^. We then substitute the asymptotic result (!25l) for 
Since the system is hnear in a^x and a^^, we can scale them such that 

axx = 1 and ayy = 7, (36) 

where 7 is the ratio ayy/ axx- We are then left with the system 



Wi[(4Wi-2)(l + 7) = 0, 
Wi7(/3-4) = 0, 



which has two solutions: 



(37) 



2 Wi-3 

^ = 0, /3 = 4Wi-2 and 7 = -^^^, = 4. (3J 



The three asymptotic results for /3 in Eqs. (153]) and (I5S|) are plotted in Figure [HI 
together with numerical results for = 10.0 and = 100.0. These were 
obtained by integrating the differential equations (ITSl) numerically from y = 
10~^° to y = 1.0 and taking the slope in a log-log representation of axx and 
ayy. Because the stress is a regular function of the extension, these exponents 
immediately carry over to the stresses around y = 0. 



13 



From the asymptotic results, we can conclude the following: although the 
stresses themselves remain finite for any Weissenberg number, for sufficiently 
large L there is a range of Wi for which j3 < 1 and stress gradients become 
infinite. This range lies roughly between Wi = 1/3 and Wi = 3/4. These 
bounds are obtained by putting /3 = 1 in Eqs. fl33l) and fl38|) . The range is 
finite, because for low Wi the velocity gradient is insufficient to cause large 
extension gradients, while for high Wi the polymers are already almost fully 
stretched well outside of the central region, which also prevents the formation 
of large extension gradients in the central region (as can be seen in Fig. [S]). 

The precise extent of the range in Wi for which the stress gradient diverges, 
depends on L^; for small it is absent entirely, as can be seen in Fig. [8] for 
= 10.0. Although the analysis in the present paper is less straightforward, 
it is fully analogous to Renardy's result for the Giesekus model [7]. 



5 Concluding remarks 

The results that we have derived here are valid for a rather artificial geometry, 
namely a purely extensional fiow on an infinite strip. However, we believe 
that the results are still relevant for more realistic fiows. The approximation 
that quantities do not depend on x is valid on the central incoming fiow 
line. If we can show that the results do not depend crucially on the exact 
boundary conditions and on the assumption of uniform extensional fiow, then 
the results that were derived above still give information on the shape of actual 
birefringent strands in fiuids that are described by the FENE-P model. 

We hope that our results for the scaling behaviour of these near-singular struc- 
tures will stimulate new experiments on the birefringence strands — we sur- 
mise that the birefringent strands seen experimentally are indeed the experi- 
mental realizations of these structures, but the data on these in the existing 
literature lack the precision to test this claim. 

Finally, we make a connection between our results and the recently observed 
instabilities in the cross-channel fiow [T] (see Fig.[T)D). There, the fully-developed 
(or being close to it) Poiseuille fiow in the inlet channels provides boundary 
conditions for the normal and shear components of the stress tensor [Tyy and 
T^y in our notation) at the infiow boundaries of the square region in the cen- 
tre of the fiow domain. These boundary conditions clearly differ from the 
special values discussed in Section 2, and thus the fiow in the central region 
is not a constant-stress solution as Eq. ([7]) and is expected to be dominated 
by strands similar to (fTTj) . Moreover, the non-zero value of the shear stress 
and the parabolic profile of the velocity at the boundaries imply that the 
actual fiow in the central domain is a combination of elongational and shear 
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components. One might then argue that the first instabihty observed in [T] cor- 
responds to switching from elongation-dominated to shear-dominated velocity 
field, while the second bifurcation would be a purely elastic instability of that 
shear-dominated flow with curved streamlines The latter are almost 

always Hopf bifurcations [23ll2l] which is consistent with the time-dependent 
flows observed in [1] at large Weissenberg numbers. Our results can be consid- 
ered to be a first step in constructing an analytic approximation to the base 
flow in the cross-channel geometry which can then be used in linear stability 
analysis. Inclusion of finite extensibility is likely to be important in studies of 
the secondary instability. 
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